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A  phase-retrieval  algorithm  is  described  and  used  to  optimize  the  phase  profile  of  a  multilevel  phase 
grating  that  achieves  a  two-dimensional  Gaussian  fan-out.  Binary-optics  technology  is  used  to  fabricate 
a  surface-relief  structure  that  implements  this  optimized  phase  profile.  The  results  of  the  fabrication 
process  are  summarized,  and  the  experimental  measurements  on  the  fan-out  element  are  presented. 
Differences  between  theory  and  experiment  arise  from  scalar-theory  approximations  used  in  the  design 
algorithm.  In  order  to  illustrate  the  scalar-theory  limitations,  a  rigorous  formulation  is  used  to  analyze 
the  diffraction  efficiency  and  the  reconstruction  error  of  a  representative  binary  surfac«-relief  phase 
grating  as  a  function  of  the  period-to-wavelength  ratio.  The  analysis  serves  to  show  qualitatively  the 
effects  of  the  period-to-wavelength  ratio  on  the  performance  of  free-space  optical  interconnects  designed 
with  a  scalar-theory  formalism. 

Key  words:  Optical  interconnects,  binary  optics,  holographic  optical  elements,  computer-generated 
holograms. 


1.  Introduction 

Free-space  optical  interconnects  are  optical  elements 
that  split  a  single  laser  beam  into  a  one-dimensional 
(1-D)  or  two-<hmensional  (2-D)  array  of  beams.'-^ 
These  interconnects  are  essential  components  in  many 
applications  of  modem  optics,  such  as  optical  parallel 
processing  and  neural  computing,®’^  in  which  it  is 
necessary  to  multiplex  light  from  a  single  source  upon 
many  detectors.  There  are  two  characteristics  of  a 
free-space  optical  interconnect;  The  degree  of  con¬ 
nectivity,  or  fan-out,  and  the  splitting  ratio.  The 
fan-out  is  a  measure  of  the  number  of  detectors 
connected  to  each  source.  Clearly,  conventional  re¬ 
fractive  optical  elements  cannot  be  used  to  achieve 
the  fan-out.  Rather,  diffractive  optical  elements, 
particularly  phase  gratings,  which  can  split  an  incom¬ 
ing  beam  into  many  beams  or  diffraction  orders,  must 
be  used.  The  splitting  ratio  describes  the  relative 
power  directed  into  each  of  the  diffraction  orders  of 
the  grating.  The  problem  at  hand  therefore  is  to 
design  a  free-space  optical  interconnect,  given  a  speci¬ 
fication  of  its  degree  of  fan-out  and  its  splitting  ratio. 
In  general,  an  exact  solution  for  the  grating  phase 
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profile  in  this  framework  does  not  exist;  however,  by 
optimizing  the  grating  phase  profile  over  one  period, 
we  can  approximate  a  specified  grating  angular  spec¬ 
trum. 

Several  phase-optimization  techniques  exist.  Sim¬ 
ulated-annealing^  and  the  Gerchberg-Saxton  phase- 
retrieval®  algorithms  are  two  of  the  most  popular. 
The  former  technique  is  well  suited  for  the  optimiza¬ 
tion  of  binary  phase  profiles,  while  the  latter  tech¬ 
nique  is  better  adapted  for  the  optimization  of  contin¬ 
uous  phase  profiles.  Here  we  employ  a  modified 
Gerchberg-Saxton  phase-retrieved  algorithm  to  opti¬ 
mize  a  multilevel  phase  profile  that  achieves  a  1-to-l  1 
Gaussian  fan-out  in  two  dimensions  (i.e.,  an  11  x  11 
euray  of  light  spots).  This  element  will  be  used  for 
an  optical  processing  system  that  serves  eis  an  artifi- 
cied  retina.’  Section  2  describes  the  algorithm  and 
presents  the  results  of  the  optimization. 

Optimized  phase  profiles  are  implemented  with 
surface-relief  structures  by  etching  away  selected 
regions  of  a  grating  substrate.  Section  3  describes 
briefly  the  binary-optics  fabrication  process  that  is 
used  to  form  such  surface-relief  structures  and 
presents  the  fabrication  results  of  the  above  multi¬ 
level  surface-relief  phase  grating.  The  measured 
etch  depths  are  considered  and  compared  with  the 
target  etch  depths  determined  from  the  phase- 
optimization  procedure.  The  experimental  measure¬ 
ments  of  the  grating  power  spectrum  are  presented  in 
Section  4.  These  measur^  values  are  compared 
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with  predicted  values  calculated  by  using  a  Fourier- 
optics  formalism  of  grating  diffraction. 

Both  the  simulated-annealing  and  phase-retrieval 
algorithms  utilize  scalar  diffraction  theory,  which  is 
an  approximate  theory  derived  from  Maxwell’s  equa¬ 
tions.  Scalar  diffraction  theory  provides  an  accurate 
description  of  grating  diffraction  for  large  period-to- 
wavelength  ratios.  Therefore  the  results  and  predic¬ 
tions  of  such  a  theory  (i.e.,  optimized  phase  profiles) 
are  consistent  with  experimentally  observed  phenom¬ 
ena  if  the  physical  period  of  the  grating  is  large 
compared  with  the  wavelength  of  the  incident  radia¬ 
tion. 

To  predict  the  diffraction  behavior  of  gratings  with 
small  period-to-wavelength  ratios,  we  must  use  a 
rigorous  formulation  of  grating  diffraction.  Of 
course,  a  rigorous  analysis  can  be  carried  out  to 
describe  all  grating  diffraction  phenomena;  however, 
such  an  analysis  is  computationally  intensive.  This 
fact  limits  the  usefulness  of  a  rigorous  theory  in  the 
design  of  free-space  optical  interconnects.  In  con¬ 
trast,  scalar  diffraction  theory  results  in  a  simple 
Fourier-transform  relation,  which  is  ideally  suited  for 
numerical  computation,  between  the  grating  transmit¬ 
tance  function  and  the  grating  angular  spectrum. 
Therefore  to  maintain  the  integrity  cf  the  solutions 
obtained  from  phase-optimization  techniques  based 
on  scaleu"  diffraction  theory,  it  is  important  to  investi¬ 
gate  the  conditions  under  which  scalar  diffraction 
theory  can  be  used  to  accurately  describe  the  diffrac¬ 
tion  characteristics  of  actual  surface-relief  phase  grat¬ 
ings. 

In  Section  5  we  apply  a  rigorous  formulation  of 
grating  diffraction®  to  analyze  the  diffraction  behav¬ 
ior  of  a  representative  binary  surface-relief  structure 
as  a  function  of  the  period-to-wavelength  ratio .  This 
structure  is  a  practical  implementation  of  a  binary 
phase  profile  that,  in  the  leu-ge  period-to-wavelength 
regime,  achieves  a  1  -to-5  uniform  fan-out.  The  grat¬ 
ing  performance,  as  a  function  of  the  period-to- 
wavelength  ratio,  is  evaluated  based  on  diffraction 
efficiency  and  reconstruction  error. 

2.  Phase-Optimization  Aigorithm 

A.  Background 

In  the  scalar  limit,  the  grating  transmittance  func¬ 
tion  and  the  grating  angular  spectrum  are  related 
through  a  Fourier  transform.  In  designing  free- 
space  optical  interconnects,  constraints  on  both  the 
transmittance  function  (i.e.,  phase-only)  emd  the  angu¬ 
lar  spectrum  (i.e.,  a  desir^  power  spectrum)  are 
imposed.  The  technique  of  phase  retrieval,®  first 
developed  by  Gerchberg  and  Saxton,  accounts  for 
these  constraints  in  an  iterative  Fourier-transform 
algorithm  for  the  optimization  of  the  grating  phase. 

The  Gerchberg-Saxton  phase-retrieval  algorithm 
results  in  a  phase  profile  that  varies  continuously  and 
is  extremely  difficult  to  fabricate.  However,  this 
continuous  phase  profile  can  be  well  approximated  by 
quantizing  a  modulo-2'n'  version  of  the  continuous 
phase  profile  into  2^  phase  levels,  where  L  is  a 


positive  integer.  This  modulo-2Tr  quantization  oper¬ 
ation  (MOD2,:/  Qm)i  where  M  denotes  the  number  of 
phase  levels  used  in  the  stepwise  approximation,  is 
the  basis  of  the  binary-optics  fabrication  process®  and 
can  be  incorporated  directly  into  the  Gerchberg- 
Saxton  phase-retrieval  algorithm.  This  modified 
phase-retrieval  algorithm  permits  optimization  of 
phase  profiles  that  can  in  practice  be  fabricated  with 
existing  technologies.  A  block  diagram  of  the  modi¬ 
fied  phase-retrieval  algorithm  is  shown  in  Fig.  1. 
A  similar  quantization  scheme  has  been  used  to  code 
computer-generated  holograms.^®  The  incorpora¬ 
tion  of  a  quantization  operation  into  the  phase- 
retrieval  algorithm  can  potentially  cause  the  algo¬ 
rithm  to  stagnate  if  the  number  of  quantization  levels 
is  small.  We  employed  a  quantization  of  16  phase 
levels  and  had  no  problems  with  the  algorithm  stag¬ 
nating. 

The  object  space  corresponds  to  the  grating,  and 
the  image  space  corresponds  to  the  grating  angular 
spectrum .  The  grating  transmittance  is  represented 
by  amplitude  distribution  t(l,  p )  and  phase  distribu¬ 
tion  i>(l  p).  Similarly,  the  grating’s  angular  spec¬ 
trum  is  described  by  amplitude  and  phase  distribu¬ 
tions  T(n,  m)  and  0(n,  m),  respectively.  The  basic 
idea  of  using  the  algorithm  is  to  Fourier  transform 
back  and  forth  between  these  domains,  while  satisfy¬ 
ing  the  constraints  in  one  domain  before  returning  to 
the  other  domain.  As  shown,  the  object  domain 
constraint  is  that  the  grating  transmittance  have  unit 
magnitude  (i.e.,  phase  only),  while  the  image  domain 
constraint  is  that  the  magnitude  of  the  angular 
spectrum  be  Gaussian.  The  Gaussian  shape  im¬ 
posed  on  the  magnitude  of  the  angular  spectrum  is 
represented  by  the  function 

G[n,  m)  -  exp -[u{n'^  +  m^)],  (1) 

where  the  parameter  <7  determines  the  desired  width 
of  the  Gaussian  fan-out  distribution. 

The  algorithm  starts  in  the  object  domain  with  an 
initial  guess  of  the  grating  transmittance  function, 
to(l,p).  Since  the  grating  is  to  be  phase  only,  the 
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Fig.  1.  Block  diagram  of  modified  phase-retrieval  algorithm. 
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estimate  is  typically  implemented  as  a  random  config¬ 
uration  of  phases, which  we  denote  <I>o(^  p).  where 
I  and  p  are  the  array  indices  that  span  a  single  period 
of  the  phase  grating.  The  initial  phase  configuration 
that  was  used  is  obteiined  by 


[/  NY  I  NY 
d>o(/,p)  =  2ir3  exp|-a  K  “  ^  I  +  1^  “  “2  ) 


where  N  is  the  number  of  array  elements  in  each 
dimension.  The  a  parameter  determines  the  initial 
duty  cycle  of  the  phase  grating,  while  p  determines 
the  initial  phase  modulation  depth  of  the  grating. 
The  initial  choices  of  a  and  p  affect  the  convergence  of 
the  algorithm  as  well  as  the  accuracy  of  the  resulting 
solution.  Various  combinations  of  initial  values  were 
tried  until  an  acceptable  solution  was  found.  In  our 
optimization  routine  the  initial  starting  values  of  a  = 
10.0  and  p  =  3.0  resulted  in  the  solution  described  in 
Subsection  2.B.  The  size  of  the  array  was  N  x  N 
with  N  =  128,  and  the  continuous  phase  profile  was 
quantized  into  M  =  16  levels.  Also,  the  algorithm 
was  terminated  after  100  iterations  ( ~  20  CPU  min. 
on  a  VAX  3600  computer).  The  quality  of  the  solu¬ 
tion  was  evaluated  based  on  two  figures  of  merit: 
the  total  diffraction  efficienQr  of  the  grating  ti<  and  the 
maximum  reconstruction  error  ARmax- 
The  total  diffraction  efficiency  represents  the  frac¬ 
tion  of  the  total  input  power  contained  within  the 
diffraction  orders  of  interest  and  is  defined  (in  one 
dimension  for  simplicity)  as 

N/2+K 

m=^2-K 

- ’  (3) 

2  Pi'n] 

m*0 

where  P[(iV/2)  -t-  q]  =  ti,  is  the  power  contained  in 
the  qth  diffraction  order,  or,  equivalently,  the  diffrac¬ 
tion  efficiency  of  the  qth  diffraction  order,  and  K 
represents  the  maximum  diffraction  order  of  interest. 
For  example,  /f  =  5  for  a  1-to-ll  fan-out  grating. 
In  general,  one  seeks  to  maximize  ti,  since  it  repre¬ 
sents  the  total  amount  of  useful  power  directed  in  the 
desired  directions. 

The  reconstruction  error  AiZ[q]  for  each  order  g, 
where  q  is  an  integer  in  the  range  [  -ff,  fif  ],  is  obtained 
(again,  in  one  dimension  for  simplicity)  by 


AR[q]  = 


where  is  the  diffraction  efficiency  of  the  qth 
diffraction  order  of  the  optimized  grating,  is  the 
ideal  diffraction  efficiency  of  the  qth  diffraction  order, 
is  the  total  diffraction  efficiency  of  the  optimized 
grating,  and  is  the  ideal  total  diffraction  effi¬ 


ciency.  As  in  Eq.  (3),  is  defined  as 

N  2+K 

-^ideali^] 

,  m=N  2-K 

=  -1^7^ - (5) 

^ideajf^  1 

m=0 

where  Pideai[^]  is  the  ideal  Gaussian  power  spectrum. 
We  can  then  take  the  maximum  reconstruction  error 

APmax  =  max  AP[q]  (6) 

Q 

among  all  the  orders  as  our  figure  of  merit.  The 
reconstruction  error  quantifies  the  degree  to  which 
the  grating  power  spectrum  approximates  the  ideal 
Gaussian  power  spectrum. 

Therefore  the  goal  in  designing  frce-space  optical 
interconnects  is  to  achieve  high  total  diffraction 
efficiencies  and  low  reconstruction  errors.  The  re¬ 
sults  of  the  phase  optimization  are  presented  below. 

B.  Results 

Figures  2(a)  and  2(b)  plot  the  ideal  2-D  Gaussian 
power  spectrum  and  the  power  spectrum  of  the 
optimized  multilevel  phase  grating,  respectively. 
The  Fig.  3  plot  slices  through  the  peaks  of  the  power 


Fig.  2.  2-D  power  spectra:  (a|  ideal  Gaussian,  (b)  optimized 
grating. 


2504  APPUED  OPTICS  /  Vol,  32,  No  14  /  10  May  1993 


spectra  in  Fig.  2  (each  normalized  to  the  respective 
total  diffraction  efficiency).  As  shown,  the  agree¬ 
ment  is  quite  good.  The  maximum  reconstruction 
error  among  these  diffraction  orders  is  0.0385,  and 
the  total  diffraction  efficiency  of  the  optimized  grat¬ 
ing  is  0.85. 

These  results  are  for  a  quantization  into  16  equally 
soaced  phase  levels.  In  general,  lower  reconstruc¬ 
tion  errors  and  higher  diffraction  efficiencies  can  he 
obtained  for  finer  quantization.  A  trade-off  exists, 
however,  if  finer  quantization  is  desired.  In  particu¬ 
lar,  the  fabrication  of  the  corresponding  surface-relief 
structure  becomes  exceedingly  difficult  as  the  degree 
of  quantization  increases.  S^ion  3  outlines  briefly 
the  grating  fabrication  process  and  presents  the 
results  of  the  fabrication  of  the  above  multilevel 
surface-relief  structure. 

3.  Grating  Fabrication 

A.  Background 

Figure  4  summarizes  the  binary-optics  fabrication 
process  used  to  fabricate  the  Gaussian  fan-out  ele¬ 
ment.  Optical  lithography  and  reactive-ion  etching, 
technologies  borrow^  from  the  semiconductor  inte¬ 
grated  circuit  industiy,  are  used  in  forming  the 
surface-relief  structure.  To  fabricate  a  surface- 
relief  structure  with  2^  levels,  we  must  design  L 
lithographic  masks,  and  L  lithography-etching  steps 
are  required,  as  shown.  For  L  >  1  an  optical 
alignment  is  performed  between  each  lithography — 
etching  step. 

Binary  surface-relief  structures  (i.e.,  L  =  1),  how¬ 
ever,  can  be  fabricated  without  optical  alignment, 
which  becomes  increasingly  difficult  as  the  degree  of 
quantization  increases  and  the  feature  sizes  decrease. 
For  this  reason  binary  gratings  are  relatively  easier  to 
fabricate  than  multilevel  gratings.  However,  as  men¬ 
tioned  earlier,  the  relative  ease  of  fabrication  comes 
at  the  expense  of  higher  reconstruction  errors  and 
lower  diflhiction  efficiencies.  A  more  detailed  expla- 


(d) 


Fig.  4.  Binary-optics  fabrication  process:  (a)  optical  lithography, 
(b)  development  of  photoresist,  (c)  reactive-ion  etching,  (d)  removal 
of  photoresist. 

nation  of  the  fabrication  process  for  both  binary  and 
multilevel  surface-relief  structures  can  be  found  else¬ 
where.® 

Two  errors  that  are  generally  present  in  the  fabrica¬ 
tion  of  both  binary  and  multilevel  surface-relief  grat¬ 
ings  are  etch-depth  and  linewidth-reproduction  errors. 
The  computer  modeling  of  the  effects  of  such  errors 
on  the  diffraction  characteristics  of  surface-relief 
gratings  is  a  current  area  of  research.  Using  a 
Fourier-optics  formalism,  one  can  show  that  for 
binary  phase  gratings  the  effect  of  an  etch-depth  error 
is  to  change  the  relative  ratio  of  the  zero-order 
diffraction  efficiency  to  the  higher-order  diffraction 
efficiencies,  while  keeping  the  relative  ratio  of  diffrac¬ 
tion  efficiencies  among  the  higher  orders  (i.e., 
where  k  ^  Q,l  ^  0,  emd  k  ^l)  constant.  For  binaiy 
gratings  that  implement  free-space  optical  intercon¬ 
nects,  as  the  degree  of  fan-out  increases,  the  sensitiv¬ 
ity  of  the  zero-order  diffraction  efficiency  also  in¬ 
creases.  For  multilevel  gratings  an  etch-depth  error 

also  primarily  affects  the  zero-order  diffraction  effi¬ 
ciency;  however,  the  relative  ratio  of  diffraction  effi¬ 
ciencies  among  the  higher  orders  does  not  in  general 
remain  constant. 

Note  that  the  above  statements  are  based  on  a 
Fourier-optics  formalism  of  grating  diffraction;  conse¬ 
quently,  they  can  be  used  reliably  to  analyze  the 
diffraction  behavior  of  actual  surface-relief  gratings 
when  the  period-to-wavelength  ratio  is  large.  The 
diffraction  behavior  of  a  representative  binary  surface- 
relief  grating  as  a  function  of  the  period-to-wave- 
length  ratio  is  considered  in  Section  5. 
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Fig.  5.  Portions  of  lithographic  masks  for  16-level  optimized  grating;  (a)  mask  1,  (b)  mask  2,  (c)  mask  3,  (d)  mask  4. 


B.  Results 

Four  lithographic  masks  were  designed  and  used  to 
fabricate  the  optimized  16-level  surface-relief  grating 
described  above.  Portions  of  these  masks  are  shown 
in  Fig.  5.  A  scanning  electron  micrograph  of  the 
fabricated  grating  is  shown  in  Fig.  6.  The  grating 


period  is  128  pm.  Table  1  summarizes  the  target 
etch  depths  and  etch  errors  for  each  etching  step  in 
the  fabrication  process.  Actual  etch  depths  were 
determined  from  averaging  measurements  taken  at 
several  locations  across  the  grating  structure.  As 
shown,  the  measured  etch  depths  are  within  2.6%  of 
the  target  values. 

4.  Experimental  Results 

A  He-Ne  (\  =  0.6328  pm)  laser  was  used  to  illumi¬ 
nate  the  grating  in  the  experimental  setup.  Figure  7 
plots  the  measured  and  expected  power  spectra  of  the 
multilevel  grating  (for  the  orders  of  interest).  The 
expected  power  spectrum  was  calculated  by  using  a 
Fourier-optics  formalism  and  accounts  for  the  etch- 
depth  errors  in  the  grating  fabrication  described 


Table  1 .  Target  Etch  Depths  and  Etch-Depth  Errors 


Target  Etch  Depth 

Etch-Depth  Error 

Mask 

(nm) 

(nm) 

Fig.  7.  2-D  power  spectra  of  a  multilevel  grating:  (a)  measured, 
(b)  expected. 

above.  Comparing  the  spectra  in  Figs.  2(b)  and  7(b), 
we  see  that  the  primary  effect  of  an  etch-depth  error 
is  to  increase  the  zero-order  diffraction  efficiency. 

Figure  8  plots  cross  sections  through  the  peaks  of 
the  2-D  power  spectra  (each  normalized  to  the  total 
diffraction  efficiency)  shown  in  Fig.  7.  With  the 
exception  of  the  zero  order,  the  two  spectra  are  in 
agreement.  Table  2  summarizes  the  reconstruction 
errors  for  the  orders  shown  in  Fig.  8. 

The  relatively  large  reconstruction  error  for  the 
zero  order  is  consistent  with  the  observation  that  the 


_  --  .4 

i/Bmcooo  croe. 

Fig.  8.  Cross  section  of  2-D  power  spectra  shown  in  Fig.  7: 
measured  (solid  curve),  expected  (dashed  curve). 


Table  2.  Reconstruction  Errors  tor  Measured  Power  Spectrum 


Q 

A/Jlgl 

-5 

0.208 

-4 

0.177 

3 

0.00230 

-2 

0.0897 

-1 

0.0427 

0 

0.178 

1 

0.0150 

2 

0.0125 

3 

0.0794 

4 

0.110 

5 

0.0724 

actual  zero-order  diffraction  efficiency  of  surface- 
relief  gratings  is  generally  higher  than  the  correspond¬ 
ing  zero-order  diffraction  efficiency  predicted  by  using 
a  Fourier-optics  formalism.  Recall  that  Fourier 
theory  (i.e.,  scalar  diffraction  theory)  provides  reli¬ 
able  estimates  of  grating  diffraction  behavior  for  large 
(i.e.,  approaching  infinity)  period-to-wavelength  ratios. 
Rigorous  analyses  of  surface-relief  gratings  have  been 
performed  that  investigate  the  deviations  of  the 
diffraction  behavior  from  that  predicted  by  a  Fourier- 
optics  analysis.12,13  Section  5  presents  the  results  of 
such  an  analysis  for  a  representative  binary  surface- 
relief  grating. 

The  reconstruction  errors  for  the  negative-fifth 
and  negative-fourth  orders  are  also  relatively  high. 
The  absolute  deviation  in  power,  however,  is  quite 
small.  For  certain  applications  such  as  optical  switch¬ 
ing,  in  which  a  certain  threshold  power  must  be 
reached,  it  is  the  absolute,  as  opposed  to  the  frac¬ 
tional,  deviation  in  power  that  is  the  more  critical; 
however,  other  applications  such  as  neural  comput¬ 
ing,  in  which  the  relative  power  contained  in  each 
order  is  tailored  to  implement  a  set  of  interconnection 
strengths,  require  small  fractional  deviations  (i.e., 
low  reconstruction  errors).  Therefore  the  informa¬ 
tion  conveyed  by  the  concept  of  the  reconstruction 
error  must  be  interpreted  with  an  application  in 
mind. 

Another  feature  of  Fig.  8  is  that  the  measured 
spectrum  has  a  slight  as3mimetry.  In  theory  the 
power  spectrum  must  be  S3nnmetric  since  the  surface- 
relief  structure  itself  is  symmetric  (e.g.,  see  litho¬ 
graphic  masks  in  Fig.  5).  Therefore  we  can  conclude 
that  the  surface-relief  structure  is  not  symmetric. 

The  occurrence  of  fabrication  errors  such  as  nonuni¬ 
form  etch  depths,  misalignment,  and  nonvertical 
sidewalls  are  all  plausible  explanations  for  the  ob¬ 
served  asymmetry  in  the  measured  power  spectrum. 
The  lack  of  computer  models  for  analyzing  the  effects 
of  such  fabrication  errors  on  the  power  spectra  of 
diffraction  gratings  along  with  the  fact  that  it  is 
extremely  difficult  to  experimentally  measure  and 
determine  the  nonuniformities  in  a  surface-relief 
profile  owing  to  such  fabrication  errors  precludes  a 
detailed  analysis  of  the  experimentally  measured 
power  spectrum. 
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An  important  point  to  note  is  that  the  minimum 
feature  size  on  the  surface-relief  structure  shown  in 
Fig.  6  is  of  the  order  of  the  wavelength  of  light. 
Recall  that  scalar  diffraction  theory  breaks  down 
when  the  dimensions  of  the  diffracting  features  of  the 
grating  approach  the  wavelength  of  light.  Therefore 
the  Fourier-optics  formalism  used  to  determine  the 
expected  power  spectrum  in  Fig.  2(b)  may  not  model 
reliably  the  diffraction  characteristics  of  the  fabri¬ 
cated  grating.  This  fact,  however,  does  not  explain 
the  observed  asymmetry  since  a  rigorous  analysis  of 
the  fabricated  grating  would  also  yield  a  sjonmetric 
power  spectrum,  owing  to  the  sjmimetry  of  the  sur¬ 
face-relief  profile. 

As  shown  above,  deviations  exist  between  the  mea¬ 
sured  and  predicted  power  spectra  (see  Figs.  2  and  8). 
These  deviations  can  be  explained,  in  part,  by  the  fact 
that  the  predicted  spectrum  was  calculated  by  using  a 
scalEU*  diffraction  theory.  A  scalar  theory,  unlike  a 
rigorous  electromagnetic  wave  theory,  does  not  ac¬ 
count  for  the  polarization  of  the  incident  light. 
When  the  grating  period  A  is  large  compared  with  the 
wavelength  of  the  incident  light  \,  the  grating  diffrac¬ 
tion  behavior  is  independent  of  the  incident  polariza¬ 
tion.  In  such  a  situation,  scalar  diffraction  theory 
can  be  used  to  obtain  reliable  estimates  of  the  grating 
diffraction  behavior.  However,  when  the  dimen¬ 
sions  of  the  diffracting  features  of  the  grating  ap¬ 
proach  the  wavelength  of  the  incident  light,  polariza¬ 
tion  effects  dominate  and  therefore  must  be  taken 
into  account. 

For  our  multilevel  surface-relief  grating  (see  Fig.  6) 
the  minimum  feature  size  is  1  nm.  Therefore  scidar 
diffraction  theoiy  cannot  be  expected  to  give  reliable 
estimates  of  the  diffraction  behavior  of  our  grating. 
For  this  reason  we  proceeded  to  investigate  the 
behavior  of  surface-relief  gratings  as  a  function  of  the 
period-to-wavelength  ratio  A/\.  Presently,  com¬ 
puter  models  that  can  analyze  the  diffraction  behav¬ 
ior  of  arbitrary  2-D  surface-relief  gratings  do  not 
exist.  However,  models  exist  that  cam  amalyze  arbi¬ 
trary  1-D  grating  structures.  For  simplicity  we  chose 
to  amalyze  a  representative  binary  surface-relief  grat¬ 
ing.  The  results  of  this  amailysis  au-e  presented  in 
Section  5. 

5.  Scalar-Theory  Design  Limitations 

Optimized  phase  profiles  for  free-space  optical  inter¬ 
connects  au-e,  in  practice,  implemented  with  surface- 
relief  structures.  Such  an  implementation  is  a  good 
one  when  the  physical  period  of  the  surface-relief 
grating  is  lau’ge  compared  with  the  wavelength  of  the 
incident  light.  In  this  section  we  amalyze  the  diffrac¬ 
tion  behavior  of  a  binary  surface-relief  grating  as  a 
function  of  the  period-to-wavelength  ratio.  This  sur¬ 
face-relief  grating  is  a  practicad  implementation  of  am 
optimized  binary  phase  grating^'*  that  achieves  a 
l-to-5  uniform  fam-out.  A  rigorous  vector  coupled- 
wave  theoiy  of  grating  diffraction,  developed  by  Mo- 
hauraun  and  Gaylord,®  is  employed.  A  computer  code 
baused  on  this  theory  is  used  to  perform  the  anadysis. 


In  our  simulations  the  incident  medium  is  air  and  the 
grating  medium  is  fused  silica  (index  of  refraction 
n  =  1.5).  Also,  the  maximum  period-to-wavelength 
ratio  considered  was  50  for  practicad  reasons  (i.e., 
necessary  computation  times). 

We  note  that  our  anadysis  is  specific  to  the  particu¬ 
lar  grating  considered.  A  similam  analysis  for  arbi¬ 
trary  grating  structures  will  in  general  give  different 
results.  It  is  not  our  intention  to  maike  general 
statements  on  grating  behavior  as  a  function  of  the 
period-to-wavelength  ratio.  In  fact,  it  is  not  clear 
whether  such  a  generalization  can  be  made.  The 
main  purpose  of  our  anadysis  is  to  illustrate  qualita¬ 
tively  the  effect  of  the  period-to-wavelength  ratio  on 
the  performance  of  surface-relief  gratings  that  imple¬ 
ment  free-space  optical  interconnects. 

The  breakdown  of  scalar  diffraction  theory  mani¬ 
fests  itself  in  grating  diffraction  with  the  onset  of 
poWization  effects.  As  the  period-to-wavelength  ra¬ 
tio  decreases,  the  grating  power  spectrum  becomes 
sensitive  to  the  polarization  of  the  incident  light. 
In  the  limit  of  a  large  period-to-wavelength  ratio  the 
grating  power  spectrum  becomes  independent  of  the 
polarization  of  the  incident  light;  consequently,  light 
can  be  treated  as  a  scalar  quantity,  as  opposed  to  a 
vector  quantity. 

Figures  9(a),  9(b),  and  9(c)  plot  the  grating  power 
spectra  for  period-to-wavelen^h  ratios  of  2,  15,  and 
50,  respectively  (dashed  curves),  and  the  puwer  spec¬ 
trum  predicted  using  a  Fourier-optics  formalism  of 
grating  diffraction  (solid  curves).  In  these  plots, 
only  the  orders  from  -2  to  +2  are  shown  and  the 
incident  polarization  is  TE .  The  key  feature  of  these 
plots  is  that  the  diffraction  behavior  deviates  signifi¬ 
cantly  from  that  predicted  by  the  Fourier  theory 
when  the  period-to-wavelength  ratio  is  small.  Con¬ 
versely,  when  the  period-to-wavelength  ratio  is  large, 
the  diffraction  behavior  of  actual  surface-relief  grat¬ 
ings  can  be  modeled  accurately  with  a  Fourier-optics 
formalism. 

Figure  10  plots  the  maximum  reconstruction  error 
versus  the  period-to-wavelength  ratio  for  both  TE 
and  TM  incident  polarizations.  The  minimum  peri- 
od-to- wavelength  ratio  considered  was  2.  This  value 
was  chosen  such  that  at  least  five  orders  were  propa¬ 
gating  in  both  the  incident  medium  and  the  grating 
medium.  There  are  three  important  features  of  this 
plot. 

First,  as  the  period-to-wavelength  ratio  increases, 
the  TE  and  TM  curves  approach  each  other.  This 
result  is  reasonable  since  we  know  that  in  the  limit  of 
an  infinite  period-to-wavelength  ratio  the  two  cases 
are  degenerate. 

Second,  at  small  and  moderate  values  of  the  period- 
to-wavelength  ratio  the  maximum  reconstruction  er¬ 
rors  for  the  TE-  and  TM-polarized  incident  waves 
deviate  significantly.  The  peirticularly  low  value  of 
the  maximum  reconstruction  error  for  the  TM- 
polarized  case  at  a  period-to-wavelength  ratio  of  9  is 
an  unexpected  result.  The  maximum  reconstruc¬ 
tion  error  at  this  point  is  0.0287,  which  is  even  lower 
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Fig.  9.  Power  spectrum  of  a  l-to-5  fan-out  grating  for  various 
period-to-wavelength  ratios  A/X  with  TE  incident  polarization: 
(a)  2.0,  (b)  15.0,  (c)  50.0.  Rigorous  calculation  is  denoted  by  the 
dashed  curve;  the  Pourier-optics  prediction  is  denoted  by  the  solid 
curve. 


Fig.  10.  Maximum  reconstruction  error  versus  period-to-wave- 
length  ratio  for  a  l-to-5  fan-out  grating  and  TE  and  TM  incident 
polarizations. 


than  the  maximum  reconstruction  error  at  a  period- 
to-wavelength  ratio  of 50, 0.0348.  The  grating  power 
spectrum  for  an  incident  TM-polarized  wave  and  a 
period-to-wavelength  ratio  of  9  is  shown  in  Fig.  11 
along  with  the  ideal  power  spectrum  predicted  from  a 
Fourier-optics  formalism.  Although  this  anomaly 
cannot  be  generalized  to  describe  the  diffraction 
characteristics  of  arbitrary  grating  structures,  fur¬ 
ther  investigation  and  research  on  grating  diffraction 
as  a  function  of  the  period-to-wavelength  ratio  may 
enable  designers  to  utilize  occurrences  such  as  this  in 
small,  compact  optical  systems. 

The  last  point  to  note  is  that  Fig.  10  contains  no 
information  on  the  actual  amount  of  power  contained 
in  each  of  the  transmitted  diffraction  orders.  For 
example,  the  plot  shows  that,  at  a  period-to-wave- 
length  ratio  of  5,  both  the  TE-  and  TM-polarized 
cases  have  a  maximum  reconstruction  error  of  approx¬ 
imately  however.  Fig.  12  shows  that  the  power 


Fig.  11.  Power  spectrum  of  a  l-to-5  uniform  fan-out  grating  for 
A/X  =  9  and  TM  incident  polarization.  Rigorous  calculation  is 
denoted  by  the  dashed  curve;  the  Fourier-optics  prediction  is 
denoted  by  the  solid  curve. 
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Fig.  12.  Power  spectrum  of  a  l-to-5  uniform  fan-out  grating  for 
.\/k  =  5  and  TE  (solid  curvei  and  TM  (dashed  curve)  incident 
polarizations. 

spectra  for  these  two  polarizations  differ  dramati¬ 
cally. 

Figure  13  plots  the  total  diffraction  of  the  l-to-5 
uniform  fan-out  grating  as  a  function  of  the  period-to- 
wavelength  ratio  for  both  TE-  and  TM-polarized 
incident  waves.  As  expected,  the  diffraction  efficiency 
is  independent  of  polarization  for  large  values  of  the 
period-to-wavelength  ratios,  and  polarization  effects 
become  apparent  as  the  period-to-wavelength  ratio 
decreases.  These  effects,  however,  are  not  as  dra¬ 
matic  as  those  for  the  reconstruction  error. 

Figure  13  also  shows  that  high  total  diffraction 
efficiencies  eu’e  possible  for  small  period-to-wave¬ 
length  ratios.  Physically,  the  increase  in  total  diffrac¬ 
tion  efficiency  for  small  period-to-wavelength  ratios  is 
caused  by  the  number  of  propagating  (i.e.,  nonzero 
time-averaged  power)  diffraction  orders  in  the  grat¬ 
ing  medium  decreasing  as  the  period-to-wavelength 
ratio  decreases.  In  theory  it  is  possible  to  achieve  a 
total  diffraction  efficiency  of  100%  (ignoring  Fresnel 


Fig.  13.  Total  diffraction  efficiency  versus  period-to-wavelength 
ratio  for  a  l-to-5  uniform  fan-out  grating  v/ith  TE  (solid  curve)  and 
TM  (dashed  curve)  incident  polarizations. 


losses)  by  decreasing  the  period-to-wavelength  ratio. 
However,  one  must  remember  that  the  reconstruc¬ 
tion  error  increases  a  great  deal  as  the  period-to- 
wavelength  ratio  decreases.  Therefore  a  trade-off 
exists  between  achieving  a  high  total  diffraction 
efficiency  and  a  low  maximum  reconstruction  error. 
The  optimum  value  of  the  period-to-wavelength  ratio 
is  determined  by  the  particular  application  at  hand 
(e.g.,  optical  switching-thresholding,  beam  combin¬ 
ing,  or  optical  implementation  of  the  synaptic  inter¬ 
connections  of  a  neural  network). 

6.  Summary 

A  modified  Gerchberg-Saxton  phase-retrieval  algo¬ 
rithm  has  been  used  to  optimize  phase  profiles  that 
can  be  fabricated  by  using  binary-optics  fabrication 
technology.  An  optimized  multilevel  phase  grating 
that  achieves  a  2-D  Gaussian  fan-out  was  described. 
The  total  theoretical  diffraction  efficiency  of  the 
grating  is  85%,  and  the  theoretical  maximum  recon¬ 
struction  error  is  3.8%.  A  surface-relief  structure 
that  implemented  the  optimized  phase  grating  was 
fabricated.  Good  etch-depth  control  was  achieved 
( <  3%  error);  however,  the  measured  grating  power 
spectrum  lacked  symmetry.  This  fact  was  attrib¬ 
uted  to  various  fabrication  errors  such  as  nonuniform 
etch  depths,  nonvertical  sidewalls,  and  misalignment 
of  the  various  mask  levels. 

Finally,  theoretical  limitations  on  the  implementa¬ 
tion  of  optimized  phase  profiles  (i.e.,  determined  from 
phase  optimization  techniques  such  as  phase  re¬ 
trieval  or  simulated  annealing)  with  surface-relief 
structures  were  discussed.  In  particular,  the  use  of 
scalar  diffraction  theory,  valid  for  large  period-to- 
wavelength  ratios,  to  describe  grating  diffraction  was 
addressed.  The  performance  of  a  binary  surface- 
relief  grating  was  analyzed  as  a  function  of  the 
period-to-wavelength  ratio.  The  performance  was 
evaluated  based  on  the  maximum  reconstruction 
error  and  the  total  diffraction  efficiency. 
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Robert  Knowlden  and  Edward  Motamedi  for  help  in 
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